Asymptotics of a Brownian ratchet for Protein Translocation 
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Abstract 

Protein translocation in cells has been modelled by Brownian ratchets. In such models, 
the protein diffuses through a nanopore. On one side of the pore, ratcheting molecules 
bind to the protein and hinder it to diffuse out of the pore. We study a Brownian ratchet 
by means of a reflected Brownian motion {Xt)t>o with a changing reflection point {Rt)t>o- 
The rate of change of Rt is 'y{Xt — Rt) and the new reflection boundary is distributed 
uniformly between and Xt- The asymptotic speed of the ratchet scales with j^^^ and 
the asymptotic variance is independent of 7. 



s 

' — ' 1 Introduction 

> 
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Brownian motion is a model of thermal fluctuations of small particles. If a particle moves 
according to such fluctuations, it is known to undergo undirected movement, i.e. Brownian 
motion is a martingale. However, the idea to use random fluctuations in order to force 
a particle in one direction is tempting and led to the paradox of the Brownian ratchet as 



introduced by Feynman et al. (1963, Chapter 46). 



Quantitative models of Brownian ratchets for biological mechanisms were introduced in 
Simon et al. (19921 and Peskin et al. (19931. While the polymerisation ratchet serves as a 



model of growth of polymers against a barrier, we focus on the translocation ratchet, a model 
for protein transport. Here, a polymer moves in and out through a nanopore by thermal 
fluctuations. On the m-side of the nanopore, molecules may bind to the polymer and bound 
sites along the polymer are forbidden to move through the nanopore; see Figure [T] for an 
illustration. 

A prominent example of a translocation ratchet was suggested on the basis of empirical 



data for Prepro-a Factor as protein and BiP as the ratcheting molecule by Matlack et al. 



( 1999 ): after translation, proteins have to be transported into the lumen of the endoplasmatic 



reticulum (ER) where they are e.g. cut and folded in order to function properly. In general, 
such ratcheting molecules can be chaperones which bind to the translocated protein in-side 
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Figure 1: (A) Illustration of the translocation ratchet. A protein moves through a nanopore 
by thermal fluctuations. On the in-side of the pore, molecules may bind to the protein which 
prevent the protein from moving out again. (B) The state of the protein at time t. The total 
length of the protein on the m-side is Xt whereas the last bound molecule is located at Rt- 



the ER lumen. The resulting translocation ratchet has been studied and compared to data 



from Matlack et al. (1999) by Liebermeister et al. (2001) and Elston (2002). In these models, 



a discrete set of sites in-side the ER lumen is either bound or unbound with BiP. The protein 
diffuses to either side of the nanopore with equal chances. Each unbound site of the protein 
may become bound by BiP and each bound BiP molecule may dissociate from the prepro-a 
Factor. Finally, when the protein is completely in-side the ER lumen, it is released. 

In mathematical models of translocation ratchets, several parameters have to be specified: 
the distance of sites along the protein which can be either bound or unbound with ratcheting 
molecules; the rate and location of association and rate of dissociation of ratcheting molecules. 
The case of a large association rate and a fixed distance of possible ratcheting sites was studied 
in Budhiraja and Fricks ( 2006 ) . (This ratchet is similar to the polymerisation ratchet of Peskin 



et al. 1993 ) Budhiraja and Fricks obtain a Law of Large Numbers and a Central Limit 



Theorem for the speed of protein translocation if the protein diffuses through a Brownian 
motion with drift. In the present paper, we study a translocation ratchet for a continuum of 
possible ratcheting sites, in the limit of small ratcheting molecules, small dissociation rates 
and long proteins. In this model, we can compute a Law of Large Numbers (Theorem [T]) and 
a Central Limit Theorem (Theorem [2]) for the speed of protein translocation. 

The paper is organised as follows: in Section [2] we introduce our model and give the main 
results, the Law of Large Numbers (Theorem [T| and the Central Limit Theorem (Theorem 
[2]). In Section [3] we interpret our results with respect to existing literature on ratchet models. 
In Sections |4j |5] and [6] we provide the three main techniques used for the proof of Theorems [T] 
and|2j see also Remark |2.2[ We conclude with a formal proof of both theorems in Section [7} 
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Figure 2: One realization of the processes {X,Tl), X = (Xt)(>o,7^ = {Rt)t>o- The Brownian 
motion X is reflected at the boundary TZ which in turn jumps at time t at rate proportional 
to Xt - Rt. 

2 Model and Results 



We study a translocation ratchet similar to the model by Liebermeister et al. (2001) and 



Elston (20021 using the following assumptions: (i) the protein moves in and out with equal 
probabilities; (ii) the protein movement is reflected at bound ratcheting molecules; (iii) the 
protein has a continuum of sites to which ratcheting molecules can bind; (iv) the dissociation 
rate of ratcheting molecules from the protein is much smaller than their binding rate to the 
protein; (v) the ratcheting molecules are inflnitely small; (vi) the protein is infinitely long. 

In mathematical terms, we consider a time-homogeneous Markov dynamics {X,TZ) = 
{Xt, Rt)t>o, starting in {Xo,Ro) = (x,0),x > 0. Here, at time t, Xt is the length of the 
protein on the m-side and Rt is the distance of the molecule closest to the boundary as 
measured from the end point of the protein that is inside the cell. The process X = (Xt)t>o 
is reflected Brownian motion, where the reflection point at time t is Rt- This reflection point 
Rt increases with t: given that the value of {X,TZ) at time t is {Xt,Rt), Rt jumps to a 
uniformly chosen value r G [Rt,Xt] at rate ^{Xt — Rt)dr for some 7 > 0. In other words, at 
rate ^{Xt — Rt) a new ratcheting molecule binds uniformly on [Rt] Xt\ which provides a new 
reflection point for X. By this dynamics, Rt < Xt for all t > 0, almost surely. We refer to 
{X,TV) as the ^-Brownian ratchet started in x. If we want to stress the dependence on 7, we 
write {X'^,TZ'^) = {X^ , R'J)t>o for the 7-Brownian ratchet. For a realization of the processes 
(A", 7^) see Figure |2j 

The rate of jumps of TZ interacts closely with the distance of X and TZ. Since TZ is non- 
decreasing and TZ < X, the process X also tends to grow. We immediately formulate our 
main results concerning the Law of Large Numbers and the Central Limit Theorem for the 
speed of X: 

Theorem 1 (Law of Large Numbers). Let X = {Xt)t>o be the ^-Brownian ratchet started in 
X > 0. Then, 

Xt t^oo „ 

— >C^ 
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almost surely, where 



r(2/3)/37y/3 



and r(.) is the gamma function. 

Theorem 2 (Central Limit Theorem). Let 7 > and X = {Xt)t>Q be the j-Brownian ratchet 



started in x >0, as in (2.1 ) and X some A^(0, 1) -distributed random variable. Then, there 



is fj > which is independent of 7 such that 



— C^t t^oq 



where ' denotes convergence in distribution. 



>X, 



Remark 2.1 (Numerical values for and o"). The numerical value for the speed of the 
7-Brownian ratchet is = W^ii^Y^^ ~ 0.4592487i/3_ 



7 ~ r(i/3) \ i J 

For the numerical value of a, note that Xt is a reflected Brownian motion for 7 = 0. Since 
a does not depend on 7 (as long as 7 > 0), it is tempting to conjecture that a = y^l — 2/7r ~ 
0.60281, the asymptotic standard deviation for reflected Brownian motion. Clearly, this 
heuristics would require an interchange of limits, t — > cxd and 7^0. While we could not 
make this interchange rigorous, simulations support this conjecture for the value of a. 

Remark 2.2 (Main steps in the proofs). The proof of Theorems [T] and [2] relies on several 
ingredients which we will provide in Sections |4j [5] and [6j respectively. 

First (see Section |4|, the process {X,TZ) can be constructed graphically by a one- 
dimensional (non-reflecting) Brownian motion and an independent Poisson process on 
[0; 00) X M with intensity 7. From this graphical construction, it becomes clear that the 
speed of Xt, i.e., the limit of must be proportional to 7^''^, if it exists. In addition, the 



graphical construction shows that Var[X(] does not depend on 7; see Remark 4.6 

Second, to compute a candidate for the asymptotic speed Xt/t, we study the Markov 
chain {Xr„ - Rr^,Rr„ - -Rr„_i,Tn - Tn-i)n=i,2,..., where To = and Ti,T2,... are the jump 
times of TZ (see Section [5]). In particular, we show that this Markov chain has a unique 
invariant distribution and (as shown in Section [7| the asymptotic speed is the ratio of expec- 
tations of Xr„ — Rt„ and r„ — r„,-i under the invariant distribution, which are computed in 



Proposition 5.8 



Third, we find a renewal structure for the process {X,TZ) where renewal points pi,p2, ... 
are given by times t where Xt = Rt and between pn and Pn+i a jump of TZ has occurred (see 
Section [6]). Using this renewal structure, we can show existence for the speed of Xt and use 
a Central Limit Theorem for cumulative processes in order to prove Theorem [2] 



3 Applications in biology and extensions 

We describe possible extensions of the 7-Brownian ratchet some of which already appeared 
in the literature. Moreover, we give biological interpretations of our findings. 



Remark 3.1 (Review of published ratchet models and extensions). Modelling protein translo- 
cation by ratcheting mechanisms has started with Simon et al. ( 1992 1 and Peskin et al. ( 1993 1 . 
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We review several features of published models for protein translocation and hint to exten- 
sions of our mathematical model. Throughout, we assume that Xt is the length of the protein 
in-side and Rt is the position of the reflection point at time t. 



Dissociation of ratcheting molecules: In the original approach of Simon et al. (19921, there 



is a finite set of equally spaced ratcheting sites along the protein. There are two rates which 
describe binding and dissociation of ratcheting molecules from the protein. In particular, the 
case of infinite rates with a constant ratio can be studied, which leads to an effective speed by 
assuming that each ratcheting site has a certain probability of being bound independent of all 
others. This is also the limiting case of fast binding and unbinding as studied in |Ambjornsson 
et al. (2005). These approaches lead to the following mathematical description: 



(i) If TZt is the set of positions of bound ratcheting molecules at time t, let the 
reflection point of Xt he Rt := maxTZt- In addition, a point r with < r < Xt is 
added at rate pdr and an existing point r G T^t is taken from TZt at a rate a. Here, 
cr describes the rate of dissociation of ratcheting molecules and the 7-Brownian 
ratchet is the special case o" = 0. 

Protein movement against a force: Usually, proteins are present in folded states inside a cell. 
Travelling through a nanopore requires unfolding of these proteins. In particular, the part of 
the protein on the out-side might still be in a folded state while the part on the in-side is 



unfolded. Hence, moving in takes place against a force, as modelled by Liebermeister et al. 



(2001); see also Budhiraja and Fricks (2006 1 . The mathematical description of such ideas 



extends the 7-Brownian ratchet as follows: 

(ii) Between reflection events, {Xt)t>o is a stochastic process with continuous 
paths, not necessarily a Brownian motion. In particular, if proteins have to unfold 
out-side the nanopore, {Xt)t>Q is best chosen as a reflected Brownian motion with 
negative drift. 

State- dependent binding rate: The binding of ratcheting molecules to the protein might de- 
pend on various parameters. In the model of Liebermeister et al. (2001) ratcheting sites 



along the protein are occupied only if they are close to the nanopore. The physical reason 
for this preferred binding close to the nanopore is an interaction with the protein forming 
the nanopore. Such a phenomenon has been described for the translocation of Prepro-a Fac- 
tor (the protein) with BiP (the ratcheting molecule) and SecGlp as the protein forming the 
nanopore (Matlack et al. , [1999 1 . A compromise between ratcheting molecules which bind 
with uniform rates along the protein and only in proximity to the nanopore was considered 



by Elston ( 2002 1 . Finally, binding of ratcheting molecules might depend on the amino acid 

We suggest the following 



(2008) 



sequence of the protein as described in Abdolvahab et al. 
extension of the 7-Brownian ratchet: 

(iii) Given Xt and Rt, let the reflection point change to r, i?^ < r < at rate 
p{d{Xt — r)) for some (j-finite measure p on R_(_. The 7-Brownian ratchet is the 
special case that p is Lebesgue measure on M_|_. 

Some special cases of the extensions (i), (ii), (iii) are straight forward to analyse. Consider 
the extension (i) as an example. As will become clear in Section |4] the average time between 
jumps of IZ in the 7-Brownian ratchet scales with 7"^/^. This implies that if cr <C 7^^^ the 
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extension (i) is approximately equal to the 7-Brownian ratchet but in other cases, the speed 
of the ratchet in (i) remains to be solved. Another example is the extension (iii) in the case 
p = 7 • (5o for some 7 > 0, i.e., binding of ratcheting molecules to the protein occur at constant 
rate directly at the nanopore. For the case of BiP as ratcheting molecule and Prepro-a as 
protein, parameters estimated in Elston (2002) indicate that this is a realistic scenario. In 
this case, jump times of 7^ = {Rt)t>o are a rate-7 Poisson process. At each jump time, we 
set Rt = Xf which is also the new reflection point of the Brownian motion. Therefore, jump 
times are renewal points of such a Brownian ratchet. Moreover, starting a reflected Brownian 



motion {Bt)t>o at and if T is an exponential time with rate 7, it can be computed (Borodin 



and Salminen 



2002 



p. 333) that Bt is exponentially distributed with rate (27)^/^. This leads 



to an asymptotic velocity of (27) 
as well. 



-1/2.. 



(7/2)1/2. ^ Central Limit Theor em can be proved 



Remark 3.2 (Biological interpretation of our results). In our model we use the asymptotics 
of an infinitely long protein and study properties of the velocity of the Brownian ratchet. 
In practice more relevant are of course a finitely long protein and the time the protein is 
completely on the m-side. For long proteins both approaches are similar: if T^ := inf{t : 
Xt > x} in our model, Tx/x ~ l/C-y with C-y from (2.1) for large x. 



Consider the case that the rate 7 for binding of ratcheting molecules to the protein, is 
proportional to the concentration a of ratcheting molecules on the in-side of the nanopore. 
As Theorem [1] shows, the speed of translocation, Xt/t, is mainly determined by 7-^/^ which 
is proportional to a^/"^. That is, to double the speed of translocation requires an eight times 
higher concentration of ratcheting molecules. In contrast, consider a ratchet model described 



at the end of Remark 3.1 where ratcheting molecules bind to the protein preferably in close 
proximity to the nanopore. We observed that the speed of the ratchet is (7/2)^/^, so the 
speed of the ratchet is proportional to o^/^. This means that the latter ratchet model uses 
the existing ratcheting molecules more efficient (i.e. protein translocation is faster) if the 
concentration a is large but the 7-Brownian ratchet is more efficient when a is small. It will 
be interesting to see if real biological systems tend to behave like the more efficient model. 

Remark 3.3 (Edwards model). One mathematical model similar in spirit to the Brownian 
ratchet is the Edwards model (see van der Hofstad et al. 1997). Let {Bt)t>o be standard 
Brownian motion and let P denote its distribution on path space. Edwards' model is a 
transformed measure on the path space defined for 7 > by 



dP^ 

Ilp 



1 



exp 



{-7 / L{T,xfdxX 



T > 0, 



where L{T,x) is the local time in x up to time T and Zj, is the normalising constant. It 
can be seen that f^L(T,x)'^ dx is the self intersection local time. Thus, the new measure 
discourages self intersections of paths. For the Brownian ratchet, such intersections are also 
discouraged by the jumping reflection boundary. 



The weak Law of Large Numbers is proven in ( Westwater 1985 ). In van der Hofstad et al. 



(1997) the Central Limit Theorem in the following form is proven: For every 7 G (0, 00) 



lim p:^, 



< c 



a' 



'2tt J- 



c 



for aU C G 



As in the case of the Brownian ratchet, the asymptotic speed is of the form C* = 6*7^/^ 
for some absolute constant b* . In addition, a* does not depend on 7 as well. The constants 
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Figure 3: One realization of the graphical construction. The shown realization of the graph- 
ical construction leads to the same path of {X,Tl) as shown in Figure [2j 

b* and a* can be expressed in terms of the largest eigenvalue of a certain Sturm-Liouville 



operator. In van der Hofstad (19981, Theorem 3, it is shown that b* G [1.104,1.124] and 
a* £ [0.60,0.66]. In the case of Brownian ratchet we have C^ = bj^^^ where b 0.459248 
and the conjectured value of cr is a ~ 0.60281. 

4 Graphical construction and some applications 

In this section we give a graphical construction of the Brownian ratchet; see also Figure |3j 
From that construction we deduce a useful scaling (Proposition 4.5 ) that will be needed in 
the sequel. 

4.1 The construction 

Definition 4.1. For 7 > let N'^ be a Poisson point process on [0; 00) x M with intensity 
measure jX'^{dt, dx), where denotes the Lebesgue measure on [0; 00) x M. Let B = {Bt)t>Q 
be an independent Brownian motion on M starting in x G M. We define tq, 5o, fi, ^i, f2, >S'2, . . . 
recursively by tq = 0, 5*0 = 0, 



and 



fn+i := inf {t > fn : ({*} x [Bt A Bt V 5„,])n7V^ ^ 0} 
I second coordinate of the unique element of 



5, 



n+l 



\{{fn+l}><[B,^^,^Sn^Bf^^,ySn]) 



(4.1) 



(4.2) 



That is, Tn+i is the first time after r„ when there is a point of the Poisson process between 
the graph of the Brownian motion and S'„,, and Sn+i is the second coordinate (the space 
coordinate) of that point. 
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For t G [Tn,fn+i) we set 

n 

Rt = ^\S^-Si^i\ and Xt = Rt + \Bt - Sn\ (4.3) 



1=1 



and {X,n) := (Xt,iii)t>o. 

Lemma 4.2. The process (J\!,TZ) is a j-Brownian ratchet started in \x\. 

Proof. First note that TZ is almost surely non-decreasing and Xt > Rt for all t > by 
construction. Consider any time t S [f„;f„+i). The rate of occurrence of f„+i is j\Bt — 
Sn\ = ^{Xt — Rt)- Moreover, Rt jumps to r uniformly chosen in [Rt]Xt\ if and only if 
l^n+i — Sn\ = \r — Rt\. By homogeneity of the Poisson process r is uniform on Xt\. 
In addition, Xt behaves like a Brownian motion, reflected at Rt- □ 



Definition 4.3. We denote that process {X ,TZ) = (Xt, Rt)t>Q constructed in Definition 4.1 
by the Brownian ratchet read off from (B,N^). To stress the dependence on 7 we also write 
{X^,n^) = {X^,R'J)t>o for this process. 

Remark 4.4. It is intuitively clear from the graphical construction that the long-time be- 
haviour of X, as given in Theorems [T] and [2] and TZ, are identical. We will prove this fact in 



Proposition 6.7 



4.2 Scaling property 

From the graphical construction we can deduce the following scaling property. 

Proposition 4.5. Let {X"',n^) = {X^',R])t>o and {X^,n^) = {Xl,R])t>Q he Brownian 
ratchets with rates 7 > and 1, respectively, both starting in x = 0. Then, 

{XlR2)t>, ^ 7-'/'(^^./3„i?;./3,),>o- (4.4) 



Proof. We use the same notation as in Definition |4.1[ By Lemma 4.2 it suffices to show ( |4.4 
for the Brownian ratchet {X ,TV) read off from B and N'^ . Setting 




^ (7-2/3t,7-i/3x). 

we find that g{N'^) = N"^. Moreover, for the space-time version of the Brownian motion 

B, i.e. B = {t,Bt)t>o, started in 0, we have g{B) := {g{t, Bt))t>o = by the Brownian 
rescaling. Applying g to each space-time point we obtain a path of a Brownian ratchet based 



on g{B) = B and g{N^) = N . In other words, (4.4) holds for {X^ , R]) and we are done. □ 

we 



Remark 4.6. 1. From the scaling property of the Brownian ratchet. Proposition 4.5 
have 

mi = „^ ,-v355!^ = y/^^S^J = y/3 u„ mi 

t—roo t t~roo t t— >oo 7^/'^t t— »oo t 
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and 



t—>00 t t—iOO t t—>oo "f'^/'^t t—>oo t 

In particular, the second limit does not depend on 7. The existence of the limits will 
be proven in Section [7] 

2. Thanks to the scaling property of the Brownian ratchet, we may choose the most conve- 
nient value of 7 = ^ in most proofs below. Afterwards (i.e. in the proofs of Theorems [l] 
and^ we use Proposition 4.5 to obtain results for general 7. 



5 The Brownian ratchet at jump times 



In this section we prove existence and uniqueness of the invariant distribution of a Markov 



chain representing the Brownian ratchet at jump times (see Definition 5.1 ). For this, we start 



in Lemma 5.4 with studying the time to the next jump and the increment of the Brownian 
ratchet during that time starting in x in expectation. We then derive properties of the 
tails of the waiting time between jumps (Lemma 5.5). Afterwards we come to existence 



(Proposition 5.6) and uniqueness (Proposition 5.7) of an invariant distribution of the Brownian 
ratchet at jump times. We conclude the section by computing the waiting time to the next 



jump and increments for a Brownian ratchet in equilibrium (Proposition 5.8). In the whole 
section we restrict ourselves to the notationally convenient case 7 = ^- Our results can easily 



be extended to general 7 using the scaling property. Proposition 4.5 



Definition 5.1. For the Brownian ratchet {X,TZ), set tq 
jump times of TZ by n, r2, .... We set (3^, W, rj) := (y. 



n, Wn, ??n)n=l,2 



and denote the sequence of 
by 



Yn=Xr„-RT„, Wn = Rrr, - Rt„-i and 7]n = Tn-Tn-l. 



(5.1) 



Observe that for any k, (Yn,Wn,r]n)n=k+i,k+2,... depends on {Yn,Wn,rjn)n=i,...,k only through 
Ifc. In particular, (y,W,rj) is a Markov chain. 

In this section we will frequently use Airy functions. We recall basic facts concerning these 
functions first. 

Remark 5.2 (Airy functions). The Airy functions Ai and Bi are two linearly independent 
solutions of the differential equations 



u"{x) — xu{x) = 



(5.2) 



with Ai{x) 



and Bi{x) 



00. We only need properties of the Airy functions on the 



non-negative real line. For further properties and explicit definitions of the Airy functions in 



terms of integrals or Bessel functions we refer the reader to Abramowitz and Stegun (1972). 



Denoting the gamma function by F we have (see Abramowitz and Stegun 
and 10.4.5) 



1972, 10.4.4 



Ai{0) 
-Ai'iO) 



1 



Bi{0) _ 

~ 32/3F(2/3) 
Bi'{0) 1 



VS 3V3r(l/3)' 



(5.3) 
(5.4) 
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The Wronskian, 

Bi'{x)Ai{x) - Ai'{x)Bi{x) (5.5) 



which is the determinant of the fundamental matrix of the differential equation (5.2), does 
not depend on x and is given by 

w := Bi'{0)Am - Ai'{0)Bi{0) = ^^,1,^,^, = -• (5.6) 

V-ii (3)1(3) vr 

(The last equality follows from Euler's reflection formula r(2;)r(l — z) = vr/sinvrz.) 

The asymptotics of the Airy functions are well known, see e.g. [Abramowitz and Steguii 
(|1972[ p. 448-449) or | Janson| ([2007| p. 161). As x ^ 00 we have 



vr 



-1/2 



Ai{x) ~ __x-i/4g-2x3/V3^ 

(Here, as usual, a{x) ~ b(x) as x — > 00 means that a{x)/b(x) ^ 1 as x ^ 00.) The integral 
of Ai is given by 

Ai{u)du=^. (5.8) 



3 



We note that 



/•X /"OO 

Gi{x) := Ai{x) / Bi{y) dy + Bi{x) / Ai{y) dy. (5.9) 

Jo Jx 

is solution of the inhomogeneous equation 

u"{x) — xu{x) = . 

vr 

5.1 Increments for single jumps of the Brownian ratchet 

Between jumps of TZ, the process X behaves like reflected Brownian motion. A jump time of 
the ratchet can be seen as a killing time of the reflected Brownian motion. Hence, we study 
a reflected Brownian motion, killed at rate 7 times its value. Recall our assumption 7 = ^• 

Definition 5.3 (Killed Brownian motion). Let {Bt)t>o denote Brownian motion started in 
X > and let 



r = inf 



{t>0: |B,|ds>c}, 



where ^ is an exponentially distributed, independent random variable with rate 1. Then, 
B := {Bt)t>o with Bt := \Bt\ for < t < r and Bt = A for t > t and some A ^ M is 
a reflected Brownian motion, killed at rate ^\B\. We denote the probability measure of the 
Brownian motion, started in x, by Fx- 

To compute the speed of the ratchet we need the expectations of r and Br- 
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Lemma 5.4. Let B he reflected Brownian motion, killed at rate ^\B\. Then, 

27rAi{x) 



E^[Br-]=X + 



3i/6r(2/3) 



and 



E^[r] = 2TT{Gi{x) +3-'^/'^Ai{x)). 



(5.10) 



(5.11) 



P roof. First we compute the Gr een function of B following the scheme outlined on p. 18-19 



Borodin and Salminen 



(2002|. Set ay := inf{t > : Bt = y}. The diffusion process B is 



m 

regular in the sense that for each x,y >0 

Fx{ay < oo) > 0. 

Such a diffusion is called transient if for some (and then for all) x,y,x ^ y 

F^ay = oo) > 0. 

As B can be killed in each interval with positive probability it is transient. Let p(.; ., .) be the 
transition density of B with respect to the speed measure, which is given by m{dx) = 2 dx. 
In the transient case the Green function, defined by 



G{x,y) :-- 



p{t;x,y)dt, 



is finite for each x,y >0 and can be computed in terms of solutions of the differential equation 
(5.2) with appropriate boundary conditions at and oo. Two linearly independent solutions 
are given by the Airy functions Ai and Bi, see Remark 5.2 The solutions of (5.2) that we 
need are obtained as follows. For x > let 

(/)(x) = Ai{x) 

tpix) = CiBi{x) + C2Ai{x). 



The function is a decreasing solution of ( |5.2[ ) and satisfies lim2,.^oo (f'ix) = 0- The constants 
Ci and C2 have to be chosen such that the function ip is increasing and (this is the requirement 
for reflecting boundary at zero) 



V''(0) = 



(5.12) 



It is easy to see, that in order to satisfy (5.12) we must have C2 = —CiBi' (0) / Ai' (0) = CiVs 
(see (5.4) for the last equality). We may take any positive number Ci in the definition of tp 
because, as one can see in (5.14), multiplication of (p and -0 by some factors does not change 
the Green function. So we set Ci = 1 and C2 = \/3- The Wronskian (which is independent 

(5.13) 



of x) is given as in (5.6) by 

w = Tp'{0)(p{0) - V'(O)0'(O) = Bi'{0)Ai{0) - Ai!{0)Bi{0) 



1 

vr 
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Now the Green function is given by 



1 w^^'il>{y)(j){x) : < y < X. 



(5.14) 



The density of the killing measure with respect to the speed measure is hy, hence it is k{y) = y 



with respect to Lebesgue measure (see Borodin and Salminen, 2002 p. 17). Furthermore the 



killing position of the Brownian motion started in x has density 

G{x,y)k{y) 



with respect to Lebesgue measure (see Borodin and Salminen, 2002 p. 14). Now we can 
compute the expected killing position starting from x: 



y^G{x,y)dy = — [(t){x) j y'^ipiy) dy + ^{x) i y"^ (f){y) dy ] . (5.15) 
We compute the integrals separately. Using the fact that Airy functions are solutions of the 



differential equation (5.2) we have 



u^Ai{u) du = uAi"{u) du = xAi' [x) — Ai[x) + C 



Thus, 



/X rx 
u^Bi{u) du= uBi'iu) = xBz{x) - Bi{x) + C. 



(x) / y V(y) dy = Ai{x) {xBi'{x) - Bi{x) + Bi{0) + V^{xAi'{x) - Ai{x) + Ai{0)) 
Jo 



and 



tpix) / y'^(p{y) dy = {Bi{x) + V2,Ai{x)) ( - xAi'{x) + Ai{x)) . 



In the last equality we have used the asymptotic relation (5.7). Substituting in (5.15) and 
simplifying we get 

1 



¥.JBr\ = -ixhix) + Aiix)iBi({)) + ^/?,Aii{))) 
w V 



where 



h{x) = {Ai{x)Bi\x) — Ai'{x)Bi(x)) = w. 



Using (5.3) we arrive at 



W 



3i/6r(2/3) 



(5.16) 
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which shows (5.10). For (5.11 ), the expected jump time of the reflection boundary if started 
from X is given by 



oo roc 



W 

2 

w 



p{t;x,y) dtm{dy) 
G{x,y) dy 

rx 

(t){x) / 'il){y)dy + tp{x) 
Jo 

{Gi{x) + 5'^/^Ai{x)) 



4>{y)dy 



(5.17) 



27r{Gi{x)+3^^^Ai{x)), 



where we have used (5.8) for the next to last equality and the definition of Gi in (5.9). □ 

In the next section we will need boundedness of second moments of r, which is a conse- 
quence of the following lemma. Note that one also could use Kac's moment formula (see e.g. 



Fitzsimmons and Pitman 1999) to compute the second moment of r directly. 



Lemma 5.5 (Bound on tail distribution of r). There exist finite positive constants c and G 
such that for any x,t >0 

Px(t > t) < Ge-^K 
Proof. We use the notation of Definition |5.3[ Clearly, 

P^(r >t)< Po(r > t), (5.18) 

so it is enough to prove the assertion in the case x = and we omit the subscript in the 

proof. Using the Brownian scaling property {Bs)s>o = (i^^^^s/t)s>o for any t > and a usual 
change of variables, we have 



\BJds = t^/^ / \BJds 



We set B = Jq\Bs \ ds and use fs to denote the density of that random variable. By the 
definition of r 



(r >t) =P(it3/2^<^) = r fB{x)e~h'^'''dx = E[e-^ 

Jo 



(5.19) 



i.e., the tail probability P(r > t) is given by the Laplace transform of B evaluated at ^t^^'^- 
We consider the Laplace transform evaluated at A > first and then set A = ht^^'^- A series 



representation of the Laplace transform of B in terms of Airy functions goes back to Kac 



( 1946 ). Instead of using the formula for the Laplace transform directly we use the asymptotic 
expansions of fjg to study the behaviour of the Laplace transform for large t. We write 



-XB] 



e-^VB{x) dx + 



-\x 



ftsix) dx 



(5.20) 
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and estimate both integrals on the right hand side for appropriately chosen r. The asymptotic 

(5.21) 



\/vr 



expansions of are (see Janson, 2007 p. 108) 

as X ^ oo. 
and 

feix) = Ae-"/'-"' {x-^ + 0{1)) as x ^ 0. 

Here, ^ is a positive constant, a = 2\a'i\^/27 ~ 0.0783 and a'^ is the largest real zero of the 
derivative of the Airy function Ai. 

Set g{x) := 22;^ log x for a; > and ^(0) := 0. Note that g is continuous and negative on 
(0, 1). We choose r > such that g{x) > —a/2 for x £ [0,r]. Then for some A > A we have 
for all x £ [0, r] 

fsix) < Ae ^x~^ = Ae < Ae 5^. 

It is easy to see that the function x ^ \x + attains its maximum in x* = (a/A)^/^. For 
A > a/r^ we have x* <r and therefore 



e'^'^feix) dx<A 



e "^"^ 2^ dx 







< Are 



— Az* — 



= Are'^"^^^^^^^/^. 



(5.22) 



Now we estimate the second integral in (5.20). From the asymptotic expansion (5.21) and 
continuity of fs it follows that there exists a finite positive constant A such that 

Mx) < 

Completion of the square in the exponent followed by a change of variables yield 

6 dx 



OO 2 



3x^ 

2 dx 



V6 
V6 



1 



-g B 

y/TT 3 JSr + A 

V3 



a2 ^ _l(3r+A-,2 

< e^— — -e 2^ > 



e 2 = e 6 ( 1 - $ 
A 



3r + A 



-Xr 



3r + A 



3r + A 



for some constant A. Here $ denotes the distribution function of the standard Gaussian 
distribution. In the next to last step we used the well known inequality 1 — ^{x) < x~^e~^ 
Combining the last two displays we obtain 



A 



e-^'^fBix) dx < A -e 

■"^^ ' - 3r + A 



-Ar 



(5.23) 



Substitution of the estimates (5.22) and (5.23) in (5.20) shows that there exists finite positive 
constants c and C such that 



me 



< Ce-^(2A) 



2/3 



t > 0. 
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Now the assertion of the lemma follows from the last equation and (5.19) 

P(r >t)= E[e 



□ 



5.2 Existence and uniqueness of the invariant distribution at jump times 

The Markov chain (3^, W, rj) describes the Brownian ratchet at jump times. Next, we show 



existence (Proposition 5.6) and uniqueness (Proposition 5.7) of the invariant distribution for 
this Markov chain. 

Proposition 5.6. Denote by (P")„gN the distribution of {y,W,ri) (introduced in Defini- 
tion 5.1) induced by a Brownian ratchet starting with (Xq, Rq) = {x, 0). For each x > there 
exists ii, • • • (^nd a distribution that is an invariant distribution for {y, W, rj) and 



-1 



^P. 



(5.24) 



k=l 



where 



■ " denotes weak convergence of probability measures. 



Proof. If the limit of a subsequence of n ^ X]fc=o exists then it must be an invariant 
distribution of {Yn, Wn, rjn)n>i- A proof of this fact in the continuous time case, that can be 



easily adapted to the discrete time case, can be found in (Liggett, 1985, p. 11). It remains 
to prove that ln~^ Yl^~^ 



P" I is a tight sequence. This follows immediately once we 

'n=l,2,... 

prove tightness of (P")n=i,2,...- To this end it is enough to show that the first moments of Yn, 
Wn and r]n are bounded uniformly in n. 



Boundedness of the first moments of r]n follows from Lemma 5.5, For the boundedness of 
the first moments of Yi , I2 j • • • and Wi , W2 , • • . , recall that R^-^ is distributed uniformly on 
[Rt^--, Xt-„-]. Hence Yn and Wn have the same distribution and it suffices to show bounded- 
ness of the first moment of l^i , 12 , . 



For this, we recall from (5.10) that 



for 



27rAi(x) 

X H — ; — < x + c 

3V6r(2/3) - 



27r^i(0) 



(5.25) 



3V6r(2/3)' 



since Ai is decreasing in [0; 00). Let J^n = (^[Xt^, R^^ : m < n) and let (C/„) be a sequence of 
iid random variables uniformly distributed on (0,1) and independent of (J\!,TZ). Using that 
sequence we have 
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and hence 



E,[Yn+l] = E,[E4(X,„^, - RrJil - Un+l)\J'n]] 
= h,[E,[{Xr,,^,-RrJ\:Fn\ 

< ^E, [Yn + c] 

= jE4y„] + Jc 



where we have used (5.25). In other words, E^[y„] < x + c for all n = 1,2, . . . which implies 



boundedness of the first moments of Yi , I2 ) • • • • D 

The following lemma shows uniqueness of the invariant distribution of the Markov chain 
(3^,W,7?). 

Proposition 5.7. Let {y,W,r]) be the Markov chain at jump times as introduced in Def- 
inition \5.1\ based on a Brownian ratchet started in x > 0. If an invariant distribution of 
{y, VV, 7]) exists, then it is unique. 

Proof. We use the graphical construction and a coupling argument for the proof. Let A^'^ be 
a Poisson process on x M with intensity measure 7 • A^. In addition, = {Bl)t>o, i = 1,2 
are two Brownian motions, started in xi < X2, such that and B^ are independent before 
T := inf{s > : B^ = B'^} and Bj = S| for t > T. Moreover, (A'^7^*) is the Brownian 
ratchet read off from (i3*, N'^), i = 1,2. Let denote the jump times of 7^*, i = 1, 2 and set 



Y' = X i - R i W = R t - R I =t'' -t"- 

n tI. t!,i "n t* -""t , In ' n ' n- 



n-l 



for i = 1,2. In order to show uniqueness of the invariant distribution of (y, W, rj) it suffices to 
show that there is an almost surely finite stopping time T' such that {Xf —X^, R\ — R\,)t>T' = 
— X"^, Rl — R'^)t>T', almost surely. 
We denote by SI the reflection point of B^ in the graphical construction, i = 1,2, i.e., given 



t G [r* ; T^^i),Sl := 5^ where , . . . , r|, ... are as in ( |4T| ) and ( |42| ) for i = 1, 2. Since 
the hitting time of B^ and B^ is almost surely finite, it is enough to construct the coupling 
using only one Brownian motion Bf starting in a; = xi = X2. So, assuming Bq = Bq, we need 
to show that T := inf{t : Rj = R^} is almost surely finite. We set Si := Sq, i = 1,2 and 
(without loss of generality) si > S2. In the case si = S2 we have T = and we are done. 
Otherwise, there are three possibilities: 

X < S2 < Sl, S2 < X < Sl, S2 < Sl < X, 

where the first and the last are symmetric. As the Brownian motion is recurrent and the 
reflection boundaries tend to jump towards the Brownian motion, it is clear that the time 
until the Brownian motion hits one of the moving boundaries is a.s. finite. Thus, it suffices 
to consider the cases 

X = S2 < Sl and S2 < x = si. (5.26) 
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Given t G [r^; t^+i), recall that SI jumps at t^^-^ to the next point in N'^ between SI and Bl, 
i = 1,2. Hence, the processes couple if and only if both components of {Bt, Sj) and {Bt, S^) 
use the same point of the Poisson process N"'. This happens if the Poisson point, say N(t, y), 
that is used in the construction satisfies 

Bt < N{t,y) < S^ < Sj or S^ < S^ < N{t,y) < Bf (5.27) 



Consider the two cases of ( 5.26| ). In the first case we consider the (almost surely finite time of 



the) first jump of Sf and in the second the (almost surely finite time of the) first jump of S}. 
Due to the symmetry of the Brownian motion and the homogeneity of the Poisson process, 
this point is below x in the first case and above x in the second case, with probability at 
least |. In other words, (5.27) is satisfied at that jump time with probability at least ^. If 



coupling did not occur, we wait an almost surely finite amount of time until either of (5.26) 
is fulfilled again. Of course it is possible that one or both boundaries jump in-between. But 
at each jump time the boundaries get closer together. Then the processes try to couple again 
upon the next jump times of either Sj or S^ with success probability larger or equal | and 
so on. Thus, the coupling time is bounded by a geometric number of finite times and hence 
is almost surely finite. □ 

5.3 Increments under the invariant distribution 



Propositions 5.6 and 5.7 guarantee existence and uniqueness of the invariant distribution of 
the Brownian ratchet at jump times. Next, we derive the increment of the Brownian ratchet 
and waiting times between jumps in equilibrium. 

Proposition 5.8. Let vr be the invariant distribution for Yi,Y2, . . . . Start the Brownian 
ratchet in vr and denote the integral with respect to the resulting distribution by Ett. Then, 

E^[Yi] = -3Ai'{0) (5.28) 

and 

E^[ti] = 6Ai(0). (5.29) 

Proof. Let Y be distributed according to vr. Let a Brownian ratchet start in Xq = Y,Rq = 0. 
Since Y is in equilibrium, we find that the distribution of a Brownian ratchet at the first jump 
time is distributed like Y. To be more precise, let = {BY)t>o be a reflected Brownian 
motion, killed at rate ^\B\, started in Y with increments independent of Y. Furthermore r 



is the killing time of as in Definition 5.3 Then 



Y = B^_U, (5.30) 

where U is independent and uniformly distributed on (0,1), since the ratchet jumps to a 
uniformly distributed point in [0;i3^]. 
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We denote by fn the density of a random variable H. The requirement (5.30) implies 



friz) = fsYui^) = frix) fB^{u)fu(-)-dudx 
Jo Jo \uJ u 

COO poo -j^ 

fvix) / fB^{u)- dudx 

10 Jz U 

roo foo 

/ G{x,u) dudx 
lo Jz 



W 
1 

W L 



/y(x) 







{z<x\ 



x)ip{u)du+ / ijj{x)4){u) du 



+ ^{z>x} / '<p{x)(i){u) du 

J z 

50 px poo poo 

fY{x)cl){x) / il){u)dudx+ I fY{x)tl){x) I (f){u)dudx 

J z J z 

pz poo 

+ / fY{x)tp{x) / (f){u)dudx 



dx 



w 



-[I^{z) + h{z) + h{z) 



We differentiate with respect to z and obtain 

fY{x)4>{x)dx, 

poo 

I'2{z) = fY{z)i^{z) I (p{u)du, 



Thus, 



I'siz) = -(j){z) J fY{x)tp{x)dx - fY{z)%lj{z) I (j){u)du. 



friz) 



1 

w 



i^{z) I fY{x)(t){x)dx + (l){z) / fYix)ip{x)dx 



(5.31) 



At this point we see that /y is strictly decreasing on [0, oo) since (p and -0 are positive. We 
differentiate this equation two more times and obtain 



mz) 



1 

W L 
1 

W L 



i;'{z) I fY{x)cP{x)dx-^iz)fY{z)(l){z) 

+ ct>'{z) r fY{x)i,{x) dx + i;{z)fY{z)ct>{z) 



(5.32) 



^'(^) / fY{x)(p{x)dx + cl)'{z) J fY{x)'ip{x)dx 
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and 



1 r 

w 

1 r 

w 



'^"(z) I fY{x)(t){x) dx - 'ip'{z)fY{z)4>{z) 

+ ct)"{z) r fyixmx) dx + (P'{z)fY{z)i;{z) 
Jo 

zip{z) I fY{x)(j){x) dx + z(p{z) / fY{x)ip{x)dx 
-fY{z){^jy{z)^{z)-<P'{z)i;iz)) 

zMz) + -fY{z)w = (zfYiz))'. 

w 



(5.33) 



Here we have used (5.31) for the next to last equality. Integrating the last equation we get 

f^{2) = zfY{z) + c (5.34) 



for some constant c. Using ip'{0) = 0, which holds by the definition of ip, and (5.32 ) we obtain 
c = /y(0) = 0. From (5.31) we know that /y is decreasing. The unique (up to a constant 
factor) decreasing non-negative solution of (5.34) with c = is given hy (p = Ai. Thus, using 
(5.8), we have 



/y(z) = 3(/>(z) = 3Ai{z),z > 0. 
The expectation of Y is given by 



(5.35) 



y7T{dy) 



poo noo 

yfY{y)dy = 3 / yAi{y) dy = 3 Af{y) dy = -3yli'(0). (5.36) 
Jo Jo 

Using (5.11) we have 

POO 

E^[ti] = E[Ey[T]] =2tt /y(x)(Gi(x) + 3~^/^Ai{x)) dx. 

Jo 



(5.37) 



For evaluating the last integral we need 



Ai'^{x) dx = xAi^{x) - {Ai'{x)f 



{Ai'm^ 



Ai'^{x) / Bi{u)dudx 



Bi(u) I Ai (x) dxdu 



Bi{u) xAf(x) - [Ai[x)f 



du 



-uAi'^{u)Bi{u) + {Ai'{u)fBi{u) du, 



Ai{x)Bi{x) I Ai{u)dudx 



Ai{u) xAi{x)Bi{x) — Ai' {x)Bi' (x) 



Ai{u) / Ai{x)Bi{x) dxdu 



dx 



uAi'^{u)Bi{u) - Ai{u)Ai'{u)Bi'{u) + Ai{u)Ai'{0)Bi'{0) du. 
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Plugging (5.35) into (5.37) and combining the last equations, we obtain, using (5.4), (5.8) 



and the fact that the Wronskian ( |5.5[ ) does not depend on x, 

E^[ti] = 6tt(^J Ai(x)Gi{x)dx + 2,-^'^ j {Ai[x)f dx 



Gvr 



Ai'{u){Ai'{u)Bi{u) - Bi'{u)Ai{u)) du 

/•oo 

+ Ai'{0)Bi'{0) / Ai{u)du + 3~^/\Ai'{0)f 



6tt(-w r Ai'{u) du + Jyli'(0)Bi'(0) - Jyli'(0)Bi'(0)) 
V Jo 3 3 / 



6^i(0). 



□ 



6 Regeneration structure 



In this section we show that the Brownian ratchet has a renewal structure and can be seen 



as a cumulative process with a remainder term (see Definition 6.2). In order to use this 
structure we will have to bound moments of several quantities, i.e. of the time between 



renewal points, (Proposition 6.4), and the increment of the Brownian ratchet at renewal 
points (Proposition 6.6). As in the last section, we fix 7 = ^. 



6.1 The Brownian ratchet as a cumulative process 

Remark 6.1 (Cumulative processes). We recall a definition of cumulative processes from 



Roginsky ( 1994). Let (T„, Vn)n>i be a sequence of bivariate iid random variables with Ti > 
almost surely. The times Ti,T2, . . . are called regeneration times. Define a renewal process 
{SMt)t>o by 



mm 



n n 

|n : ^ Tj > t| and 'S'n, = ^ Vi. 



i=l 



i=l 



The process {SMt)t>o is called a type A cumulative process. Cumulative processes are well 



studied (see e.g. Smith 1955 Roginsky, 1994 and references therein). Most importantly, it 
is known that finite second moments of Ti and Vi are sufficient for the strong Law of Large 
Numbers and the Central Limit Theorem for {SMt)t>o- 

Definition 6.2 (The Brownian ratchet as a cumulative process). Given a Brownian ratchet 
{X,TZ) = {Xt, Rt)t>o with (Xq,Rq) = {x,0), x > 0, we define a sequence of regeneration 
times using times when Xt = Rt as well as jump times of TZ. We set 



po := inf{t >0:Xt = Rt}, 
po := inf{t > pq: Rt^ Rt-}, 



(6.1) 



and for n = 1, 2, . . . , using Rt- = \\ms^t,s<t Rs 



Pn ■■= inf{t > p„-i : Xt = Rt}, 
Pn := inf{t > Pn,: Rt^ Rt-} 



(6.2) 
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such that po < po < pi < pi < ■ ■ ■ ■ Note that if x = 0, then po = and Xp^ = 0. It is clear 
that {pn — Pn~i,Xp^^ — Xp^_-^)n>i is a sequence of bivariate iid random variables. We define 

Mf := min{n : pn > t}, 

n 

X 



i=l 

At := Xp„ + Xt 



PMt 



Then we have 

that is, A" is a type A cumulative process with remainder At. 



(6.3) 



Remark 6.3 (The regenerative structure in the proof of Theorem [T] and |2]) . As mentioned in 
Remark 16. 11 finite second moments of pi — po and Xp^ — Xp^^ are sufficient for the strong Law 



of Large Numbers and the Central Limit Theorem for Smi ■ To see that the same holds for 
X, showing Theorems [T] and |2] we shall show in Proposition 6.7 that the remainder divided 
by \/t converges to almost surely as t — > oo. 

6.2 Tail estimates and moment bounds 

In the rest of the section, we denote the law of a Brownian ratchet {X,TZ), started in {x,0) 
by ¥x (with expectation operator E^;). We use the notation of Definition 6.2 We can use 
the Law of Large Numbers and the Central Limit Theorem for cumulative processes given 
pi - po and Xp^ 



XpQ have a finite second moment. These are the main results of this section 



(Proposition 6.4 and 6.6). 



Proposition 6.4 (Moment bounds for pi — po). Let po,pi be as in (6.1) and (6.2). Then, 

Ex[{pi-pof] < oo. 



(6.4) 



The proof is based on Lemma 5.5 and the following result. 



Lemma 6.5. Let a be the stopping time when the Brownian ratchet hits the moving reflection 
boundary for the first time, i.e. 



a := inf{t > : Xt = Rt}. 
Then, there exists a finite positive constant C such that 

sup [0-2] < C. 



(6.5) 



Proof of Proposition 6^. Clearly, pi — po = {po — Po) + {pi — Po)- Since po — Po = t with r 
from Definition 5.3, finiteness of the second moment of the first summand follows immediately 
from Lemma 5.5 Moreover, given Xp^ — Rp^ = z, we have that pi — po is distributed as a for 



a Brownian ratchet started in z with a from (6.5). Hence, finiteness of the second moment 
of the second summand follows from Lemma 6.5 Thus, the proposition is proved. □ 
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Proof of Lemma 6^. Recall the jump times ri, T2, . . . oi TZ. We consider the process Z := 
{Zt^a)t>o for Zf := Xt — Rt which is started in x, and has discontinuities at ti, T2, . . . . Note 
that a is the time at which Zt hits for the first time. 

Note that Z is an autonomous Markov process: locally it behaves like Brownian motion, 
reflected at and each time there is a point in a Poisson point process on [0; oo) x M in [0, Zt\ 
the process Zt restarts from this Poisson point, i.e. at discontinuity points the process Zt 
jumps down. 

Our goal is to bound a from above by a simpler random variable a. To this end we couple 
Z io a process Z = {Zt)t>o which starts in x and has the following dynamics: the local 
increments of Z follow exactly those of Z^ i.e. Z and Z use the same underlying Brownian 
motion. At jump times r„ of Z with Zr„ < 1, we set Z^-^ = 1. In addition, Z jumps at rate 
(1 — Zt)^ to 1. At jump times r„ with Zr^ > 1, the process Z does not jump. We denote the 
first hitting time of ^ of by a. By this coupling, we have achieved the following: 

1. At joint jump times r of ^ and Z we have Z^ < Z^- = 1. Jump times r exclusive to Z 
must satisfy Zr < \ = Z-y. For a jump time r exclusive to Z^ and if Z-r- < Zr-^ we 
have Z^ < Z-j— < Z-,— = Z^-. In particular, for all t > we have Zt < Zt- 

2. The process Z jumps at rate 1. If r is a jump time of Z, then Zt = 1. 

In particular, o" < 5, almost surely, by 1., and it suffices to show that a has a finite second 
moment. Set tq = and denote the jump times of Z by Ti,f2, .... In addition, let {Bt)t>o 
be Brownian motion starting in 1, set 



and 



r = inf{t >0: Bt = 0} 



fi+i -fi,i> 0. 



It is clear that Tq, Ti, . . . are iid random variables with rate one independent of Zt. Indepen- 
dent of initial position Zq the process Z restarts from 1 at time Ti. Define p := Pi(r < Ti) 
and let M be geometrically distributed with parameter p independent of Ti,T2, . . . , i.e. 
P(M = n) = (1 — p)p"', n > 0. Note that M is the number of fail trials for the process 
Z to reach 0, started in 1, before the next jump. It follows that 

M 
1=0 

where r is distributed as r, conditioned to be smaller than Ti and hence 



IEi[t2] 



Ei[r2,r<ri] ^ Ei[T2] 2 



Mr < Ti) 



< 



p 



p 



Combining these results, 

lExk^] <IEx[5^] <IE 



M 
i=0 



M 



< E 



(2(M + 1) + (M + 1)M + -(3 + M) 
V r> 



=: C < oo 

The right hand side is independent of x, and we are done. 



□ 
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Proposition 6.6 (Moment bounds for Xp-^ — Xp^). We have 

E,[X2J < CX), 
E,[{Xp,-Xp,f]<oo 



and for each t >0 



E,[iXt - Xp,,f] < 2{E,[{Xp, - Xp,f]+E,[Xl]). 



(6.6) 
(6.7) 



(6. 



Proof. Note that \X 



Pol 



-Bpol, which already shows that 



by the second Wald identity where we have used that ^^^[po] < oo by Lemma 6.5 
Now, using the strong Markov property at time po we have 

E,[iXp, - Xp,f] = Eo[x2j = Eo[Bl] = Eo[pi] = E^ - po] < oo, 

where we have again used Wald's second identity that is apphcable by Proposition 6.4 This 
proves (6.7). 



It remains to show (6.8). We have 

E,[iXt - Xp,,^f] = E,[{Xt - j|,<,„|] + E,[{Xt - j|,>,„j 

The first summand on the right hand side can be estimated as follows 
E^iXt - Xp,)H{t^p^^y] < E,[X^l{,^p^^y] + E,[X2J 



< E,[Xlp^^] + E,[X2J < 2E,[X2J, 



(6.9) 



(6.10) 



where for the last inequality we used the fact that the stopped process (^Xf/^p^^ is a sub- 
martingale (as usual a Ab = min{a, 6}). Using a similar argument we have for the second 
summand on the right hand side of (6.9) 



E,[iXt-Xp,jH 



EoiiXt - Xp, fl{t^p^y] < Eo[X^l{t<p,}]+MX^p^] 
Eo[Xf^p,] + Eo[x2j < 2Eo[x2j = 2E,[(X,, - Xp,) 



Combining (6.9), (6.10) and (6.11) we obtain (6.8). 

Proposition 6.7 (Asymptotics of At and Xt — Rt)- We have 

Xt — Rt 



(6.11) 



□ 



At_ 

Vi 



and 



Vi 



a.s. as t — > oo. 



Proof The proof of the first assertion is an adaptation of the proof of Lemma 8 in (Smith 



1955 p. 26). In view of (6.6) it is enough to prove the first assertion for Xq = 0. Then we 



need to show that 
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We set 



Yn= sup \Xt-XpJ. 

t£lpn-l,Pn] 



Then Yi , 1^2 > • • • is a sequence of independent identically distributed random variables with 
finite second moments. To prove this we note that 



^1 < Xp^ + SUp|XjApi|. 

t>0 



By Doob's maximal inequality 



supIXmpJ' <4Eo[X;J. 

L t>0 -I 



Thus, finiteness of second moments of Yi follows from (6.7). Now 



Note that by Theorem 1 in |Doob| ( |1948| ) we have 

Mt 1 



lim — TT-,r 
t^oo t IE[pi — Po\ 



> 0, a.s. 



(6.13) 



(6.14) 



Using Lemma 7 in (Smith 1955 p. 26) we obtain Y^/^ / Mt — > almost surely and therefore 
( [6l2| ) holds. 

For the second assertion we have 



t ~ t 



a.s. as i ^ oo. 



For t G [pn-i,Pn], n > 1 we have 



0<Xt-Rt< Y„,. 



Thus, 



{Xt - Rt)l{t>po} 



t 



a.s. as t ^ oo 



follows by the same argument as before. 



□ 



7 Proofs of Theorem [T] and [2] 



In the case 7 = 0, Theorem 1 is true since Xt = \Bt\ for a Brownian motion, started in x. 
Hence, we can assume 7 > in the rest of the proof. 

We make use of the regeneration structure, set out in Definition |6.2[ We set 



^x[pi - Po], P = ^x[Xp^ - XpJ and /J^ = Var^ 



Xp^ - X, 



ipi - Po)P' 



po 



(7.1) 



Here, r, p and f3 are independent of x due to the regeneration structure. By Propositions 6.4 



and 6.6 the time and space increments pi — po and Xp-^ — have finite second moments. 
As in (6.14) we have lim^^oo ^ — r s-i^ost surely. 
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From Proposition 6.7 (6.14) and Law of Large Numbers applied to Sn it follows that 



Xt 
t 



At MtSMt 
t t Mt 



a.s. as t ^ oo. 



We will compute ^/r using the ratio limit theorem for Harris recurrent Markov chains (see 



e.g. Revuz 1984 1 . Let vr denote the invariant distribution for {y, W, ry) determined by Propo- 
sitions 15.61 and 15.71 To show that this Markov chain is Harris recurrent we have to show that 



each Borel subset of M>q, say B, having positive Lebesgue measure is visited infinitely often 
with probability 1. We cannot apply the Borel-Cantelli lemma directly, because the events 
{{Yn,Wn,r]n) € B}n>i are not independent. But events lying between different regeneration 
points are independent. For n > we set 

B„ = {3m > 1 S.th. Pn-l < Tm-l,Tm < Rn) and {Ym, Wm, Tim) G B} . 

The events Bn are independent and have the same positive probability. Thus, applying the 
Borel-Cantelli lemma to the sequence {Bn)n>o we obtain 

r,{iY^,Wn,Vn)€^BLo.) >P,(5„i.O.) =1. 

Now using the ratio limit theorem we obtain that 



r 



lim — - 



T Rt 
um — 

i^oo t 



Rt. 

um — 



lim „ — 



a.s. n 



oo. 



where for the second equality we have used Proposition 6.7 We recall that E,r[l^i] = ^^^^[yi]. 
Furthermore, recalling that we have used 7 = ^ in Section [sj and using (5.3), (5.4) as well as 
Proposition 5.8 we obtain 



r 



1/3 3M0) ^7^3V3r(2^ 
^ ^' 6Ai{0) 22/3 r(l/3) ■ 



by the scaling property. Proposition 4.5 This proves Theorem [T 
Using the CLT for cumulative processes (see e.g. Smith 1955J 
sition [6^71 we obtain that for all j; G M 



Roginsky 1994) and Propo- 



lim 



< X 



lim ] 

t^oo 



lim ] 



At 



H — — , , , < X 



/3{t/ry/^ /3{t/ry/^ 
SMt - if^/r 



< X 



$(x). 



/3(t/r)i/2 

where $ denotes the distribution function of the standard normal distribution, and /^^ and r 



are as defined in (7.1 ). Hence, Theorem 2 holds for a = (ij \/r. 
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